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■ Abstract 
^ I 

I The coupled non-radiative excitations of the electromagnetic field and pha- 

00 ■ sons in films with a quasi one-dimensional charge density wave (CDW) 

. are evaluated for P-polarization and CDW conducting axis inside the film. 

I The prominent features are two surface phason-polariton branches extend- 

ing from the CDW pinning frequency to the frequency of the longitudinal 
optical phason. These surface phason-polariton states are confined to a 
finite band of longitudinal wave numbers. Besides surface polaritons, in- 
I finite series of guided wave modes are found which extend to large wave 

^ ' numbers. These differences to usual phonon-polaritons are caused by the 

Q . extreme anisotropy of the electric CDW reponse. At finite temperatures, 

quasi-particles are thermally excited. Their dissipation leads to polariton 
damping. Significant level shifts of surface phason-polaritons and guided 
wave modes are also found. They are due to thermal dressing of the longi- 
I tudinal optical phason via quasi-particles. The zero and finite temperature 

results including the case of neglected retardation are displayed in detail. 



Keywords: A. thin films, D. charge density waves, D. optical properties. 



1 Introduction 

The microwave, far infrared, and optical properties of charge density wave (CDW) 
in quasi one-dimensional conductors (for reviews of CDW physics see [1-3]) are 
characterized by the large polarization which arises when the CDW is deformed 
along the chain direction. 

At low temperatures, the optical dielectric function for the chain direction is 
essentially that of a polar crystal - albeit with very different frequency scales due 
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to the heavy Prohhch mass [4,5]. Together with the extreme anisotropy this leads 
to pecuhar far infrared responses. 

The elementary excitations of CDW which couple directly to the electric field are 
the phasons. Their theoretical discussion has a long history [5-15]. Phasons give 
rise to phason-polaritons [15] when coupled to the electromagnetic field. 
In this paper we study theoretically surface phason-polaritons in CDW films 
which can now be fabricated [16]. 

Surface polaritons are well known in conjunction with optical phonons, plasmons 
and magnons and have been reviewed in [17-20]. The present paper is directly 
related to the work of Khewer and Puchs [21] on surface phonon-polaritons in 
isotropic films. We generalize their work to anisotropic surface phason-polaritons 
involving phason dispersion functions calculated in [15]. Uniaxial half-space po- 
lariton problems have been considered in [22-24] (cf. also [18]) for other polariton 
mechanisms and mainly for non-retarded interaction. 

We will also investigate the effects of finite temperatures when thermally excited 
quasi-particles (qp) cause surface phason-polariton damping and level shifts. 
It is excepted that surface phason-polaritons can be investigated experimentally 
by attenuated total reflection (ATR) [25], possibly also by electron energy loss 
spectroscopy (EELS) [26], and by low energy electron diffraction (LEED) [27]. 
The guided waves which we find along with the surface phason-polaritons can be 
measured in microwave cavity experiments similar to those on superconducting 
films [28]. 

2 Basic Theory 

The bulk results for the dielectric tensor e.y(q, c<j) of a CDW given as equations 
(41-44) in [15] refer to long wave lengths q <^ ^q"^ = A/vp (A: half gap, vp- 
Fermi velocity) when q points into chain direction. The amplitude coherence 
length ^0 of CDW is rather small away from T^. Typical values are ~ 3nm 
due to the large value of the half gap A. The bulk dielectric functions can thus 
be used for wave lengths q < 10® cm~^ and for film thicknesses 2d > 30nm. The 
geometry under study is shown in Fig. 1. 

The film is embedded in a dielectric medium with dielectric constant eoo- The 
chain direction of the CDW system hes in the film and is chosen as the x-axis. 

The corresponding wave number is henceforth called q. For c? — oo , the limiting 
case of a half space is obtained. We consider non-radiative electric modes (P- 
polarization) . The electric field inside the film then has components and Ey 
and the magnetic field is transverse, lying perpendicular to the conducting axis 
(B^ 7^ 0). The three-dimensional wave vector then consists of the longitudinal 
component qx = q and the transverse component qy = k = —ia. As shown in 
[15], the generalization of the solvability relation in [21] to the present anisotropic 
case is 
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Figure 1: CDW film of thickness 2d and with dielectric function ejj(q, a;) embed- 
ded in a medium of dielectric constant eoo- The conducting CDW-axis is chosen 
as X-axis. 



- ely^ + a^Syy - q^Sxx + ^iqae^y = 0. 



The matching condition for the electromagnetic fields takes on the form 



^00 



coth{ad) 
tanh(Q;(i) 



^xx^yy ^xy 



(1) 



(2) 



^^yy ~^ '^Q^xy 

The coth-equation gives the high frequency and the tanh-equation the the low 
frequency branches. In comparison to [15], we changed the notation of the di- 
electric functions: en — > Cxx, ^tm ~^ ^yy, ^it ~^ ^xy Foi' surface polaritons, a is 
positive and defines the decay length in y-direction. The oscillating guided wave 
modes have imaginary a — ik. The light line of the embedding medium is ckq = 
where ao is defined by 



an 



-eoo- 



Eliminating Cyy from (2), using (1) gives the simpler formula 

Q€-XX 
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(3) 



(4) 



for the matching condition. 
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Equations (1), (3), and (4) determine the desired dispersion relation uj — uj{q) 
by eliminating a. This is clearly a complicated task. Enormous simplifications 
occur in the limit T ^ when all qp efi^ects freeze out in the CDW dynamics. 



3 CDW Dielectric Tensor at Zero Temperature 

The relevant components of the dielectric tensor at zero temperature according 
to [15] have the simple form 



6..(g,fc;a;) = 6o + 6Ary,(g)-^^|^^^, (5) 

o2 2 2 22/\|22 

Here, 



is the dielectric constant from virtual transitions across the gap 2Aq [5]. 
In (6), K is the Thomas-Fermi wave number and 2Ao the zero temperature gap. 
ca is much larger than the lattice dielectric constant eo in chain direction. Equa- 
tion (6) is a valid form for a; < Aq. rj and rji are dispersive corrections which 
become unity for g — > [15]. The phason velocity q at zero temperature and 
for the chain direction is given in terms of Fermi velocity vp, Peierls phonon 
frequency luq, electron phonon coupling constant A, and the half gap Aq as [5] : 

Ci = ^^Vp < (7) 

The characteristic frequencies of the zero temperature theory are the frequency 
culo of the longitudinal optical (LO) phason 

^lo = ^ = |Aa;J, (8) 
and the pinning frequency cuq which usually satisfies 

(^0 -C cvlo- (9) 

This is consistent with the clean limit assumption underlying [15]. The frequency 
lulo is of order 10^^ and thus much smaller than usual optical phonon frequen- 
cies. For this reason, the transverse lattice dielectric function €t will be taken as 
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constant in the considered frequency range. This is a strong assumption in view 
of the possible appearance of bound collective mode resonances in the FIR [29]. 
A useful dimensionless measure of the wave number q for phason-polaritons then 
is (c: light velocity) 

Z = (10) 

Z is of order unity for the present purposes, i.e., q < 10^ cm~^ . From cfg^ = 
etUJ^QZ {ci / cY it is seen that all dispersive effects inside the CDW dielectric func- 
tions are irrelevant. Thus one obtains at zero temperature 



Equation (11) is identical with the dielectric function of a polar crystal, cjq playing 
the role of the transverse optical frequency uoto and that of the dielectric 
constant eoo- However, (11) refers only to the chain direction. The Lyddane- 
Sachs-Teller relation 
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LO 



e* = eA (12) 

^0 

docs not give the static dielectric constant of CDW -which does not exist- but the 
"plateau" dielectric constant e* [30] as explained in [31]. The decisive difference 
to the usual phonon polariton case is the strong anisotropy e^x 7^ and the 
different frequency scale. With these simplification the solvability relation (1) 
reduces to 



(13) 




The matching condition becomes 

^00 J coth(a(i) | e 
Q/Q I tanh(a(i) J 



XX 



a 



(14) 



Equations (3), (13), and (14) completely determine the excitation spectrum. 



4 Spectra at Zero Temperature 

We will display the phason-polariton spectra in the W — Z plane. Z is given by 
(10) and the reduced squared frequency W is defined by 

W^^. (15) 



LO 
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The existence of surface phason-polaritons requires 



e* > eoo, (16) 

which can be fulfilled. The W — Z quadrant is dcvided into several regions by 
the light lines W = etZ /cqq {a^ = 0) and W = Z {a^ = 0) plus horizontal 
lines W = 1 and W = Wq where is also zero due to vanishing of e^x- The 
quantity Wq means lUq/luIq. Fig. 2 shows this partitioning of the W — Z plane 
for €t — 2eoo- The shaded areas have > 0. They -potentially- could contain 
surface phason-polaritons. We call them region 1, 2, and 3, respectively. The 
region to the left of the light line ao = belongs to ciq < , i.e., to radiative 
solutions which we do not consider. The remaining regions I, II, and III have 
< 0. They -possibly- contain guided wave solutions. The answers to these 
questions require detailed investigations of (14), paying special attention to the 
sign of a and of exx- 

We introduce dimensionless parameters according to 

h = — , go = — , Do = ^JT^dixiLolc. (17) 

^00 ^00 

For a film of thickness 2d = 1000 nm [16], a typical value of Dq would be 0.2. 



4.1 Surface Phason— Polaritons 

Surface phason-polaritons exist only in region 1 where (14) translates into 



{yV-Wo){W-Z) \ cothX \ {\-W){W-Z) ^ ^^^^ 



\go{hZ -W){\-W) \ tanhX, "V W -W( 

For each thickness Dq? there is one surface phason-polariton branch from the 
coth-equation and one from the tanh-equation separated by the limiting curve 
for L>o ^ oo which corresponds to the half space problem. The coth-branch ends 
at the fine W — 1 and the tanh-branch at the line W — Z , ats shown in Fig. 
3. Both branches start at {Z = Wo/h, W — Wo) and have W — hZ as tangent 
there. The limiting curve for Dq — > oo is given implicitly by 

The quantity go can be much larger than unity. Then the coth-branches (one for 
each Do ) are squeezed to the ckq = line and thus vanish. 
There are no solutions in regions 2 and 3. 
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Figure 2: Partitioning of reduced frequency versus wave number plane for the 
polaritonic response of a CDW film at zero temperature. W — uj^/ujj^Q measures 
squared frequency and Z = c^q^ / {ulpet) squared wave number in chain direction. 
W = Z and W = 2Z are the light lines for CDW film and embedding medium, 
respectively. Regions 1, 2, and 3 sustain surface phason-polaritons (a^ > 0^ while 
regions I, II, and III allow for guided waves {o? < 0). The region to the left of 
W — 2Z contains the radiative solutions. Wq measures CDW-pinning. 



4.2 Guided Waves 

In the regions I, II, and III, is negative, i.e., the transverse wave number k is 
real and the modes oscillate across the film. For imaginary a, we find in region 
I a series of solutions ( p= 0,1,2,3,...) 



^ (Z -W)(l-W) TT 

Do\ „^ = p- + arctan , 



[Z-W){W-Wo) ^20) 



W-Wo ^2 \\go{hZ -W){l-W)' 

Modes with p = 2m (m= 0,1,2,...) are from the former tanh-equation and those 
with p = 2m + 1 are from the coth-equation. The p — mode is the continuation 
of the corresponding surface phason-polariton. The p > modes start linearly 
at {Z = Wo/hjW = Wq). All modes go asymptotically towards for large 
q. Fig. 4 shows these guided wave modes. The slopes of the p > modes near 
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Figure 3: Surface phason-polariton branches of a CDW film for several values of 
the scaled thickness Dq. The limiting curve L>o — > oo (dash-dotted line) separates 
the region of even solutions (to its right) from the region of odd solutions. The 
quantities go and h are the high frequency CDW dielectric constants along and 
transverse to the chains in relation to the embedding medium dielectric constant, 
respectively. 



Z = Wo/h are given by - W^o)/(-Dg(l - Wo) +p'^{7r/2f) and decrease with 
increasing order p. 

Guided waves also exist in region II but not in III. Clearly, there must be a 
continuation of the coth-surface-polariton branch of region 1 into region II. This 
is the n—1 solution of 



^ {W-Z){W-1) TT 

— jJo\ 777 777 — ~^7r + arctan , 



(W-Z){W-Wo) ^21) 



W-Wo 2 \\go{hZ -W){W -ly 

The modes with n — 2m {m = 1,2,3,.. .) are from the former tanh-equation, the 
others (n = 2ni + l. ni = 0,1,2,3,...) are from the coth-equation. All modes 
with n > 2 start at discrete points {Zm, Wm = hZ„i) on the «o = line and go 
asymptotically to the W = Z line. The points Z^ solve the equation 

/ TT ^2 



ih-l)Zm{hZ^-l)=n^[^-^J {hZ^-Wo). (22) 
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Figure 4: Guided wave modes of region I in a CDW film for a fixed (scaled) 
thickness Dq — 1 and for several mode numbers p. The p — Q mode is the 
continuation of the surface phason-polariton in region 1 (dash-dotted line). The 
quantities go and h are the high frequency CDW dielectric constants along and 
transverse to the chains in relation to the embedding medium dielectric constant, 
respectively. 



These modes, therefore, start at rather high frequencies, the lowest one at Z2 = 3 
for our choice Dq = 1 and h = 2. 

Fig. 5 summarizes all our results for the excitation spectra in CDW films at zero 
temperature. 



4.3 Neglect of Retardation 

It is interesting to consider the limit c — > 00, i.e., when retardation can be 
neglected. In this limit, only region I survives leaving only guided wave modes. 
They follow from 

l~ 1-W TT I W-Wo , , 

rW^r^2^^'''^''i g,h{X-Wr ^^=0'1'2>3,.... (23) 
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Figure 5: Survey of all polaritonic excitations of a CDW film: Surface phason 
polaritons in region 1 (full line) and guided wave modes in regions I and II (bro- 
ken lines, numbered by integers p) for a fixed (scaled) thickness Dq — 1. The 
quantities go and h are the high frequency CDW dielectric constants along (e^) 
and transverse to the chains (et) in relation to the embedding medium dielectric 
constant eoo, respectively. The limiting curve for Dq ^ oo (dotted line) is the 
surface phason-polariton of a CDW half space. 



Here, Z combines wave number q and film thickness 2d according to 

Z = '-^{qdf. (24) 

The p = curve for small Z was calculated in [15] and has the characteristic 
behaviour oo ~ ^/q of a two-dimensional plasmon for lvq = 0- The modes with 
p > 1 are low lying and are acoustic for a;o = as shown in Fig. 6. 

5 Spectra at Finite Temperatures 

At finite temperatures, the expressions for the components of the dielectric tensor 
become very complicated and contain additional frequency and wave number 
dependencies. They also depend on elastic and inelastic scattering rates of the 
qp. 
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Figure 6: Non-retarded limit of polaritonic response of CDW films at zero tem- 
perature. The wave number scale is Z — e/^q^d^/et (2d: film thickness). The 
quantities go and h are the high frequency CDW dielectric constants along (e/^) 
and transverse to the chains (et) in relation to the embedding medium dielectric 
constant eoo, respectively. The limiting curve for Dq ^ oo (dotted line) is the 
surface phason-polariton of a CDW half space. 



At sufficiently high temperatures of about Tc/2 or higher, the inelastic scattering 
becomes relevant and the corresponding formulae for e^j arc given as eq. (62) in 
[15]. Inserting these expressions into the solvability condition (1), a surprising 
number of cancellations occur. A further significant simplification is the neglect of 
dispersion in CDW functions by utilizing the argument of Sec. 3. In the end, this 
simplification amounts to the following approximations for the dielectric tensor 



eA(T) 



-yy 



1 + 
1 + 



CO n^{T) 



(25) 



Here, uJi^tiT) are the dielectric relaxation frequencies of qp for longitudinal and 
transverse directions, respectively: 
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(26) 



eA{T) ■ ■ et 

The dc-conductivities a^^\T) are activated and account for the main temper- 
ature dependence. They are exphcitly given by equation (57) in [15]. We will 
further assume a quasi-one-dimensional situation when ujt -C uii holds. Then the 
solvability relation (1) reduces to 



a 



6a(T) 



1 + 



iuJi{T) 



LO 



(T) 



fi2 



2 ^ 



(27) 



This expression must be inserted into (14) to get the spectra uo = uj{q,T). It 
is noted that uj becomes complex at finite temperatures and that the imaginary 
part of UJ measures polariton damping. 

Temperature dependencies appear not only in cui but also in culq,^'^, ^-nd ea- 
The temperature dependent LO phason frequency is given in terms of its zero 
temperature value culq and the backflow parameter bo as 



UJ 



LO 



(T) 



cKTW 



The expression for fl"^ is 



(- 



An 



6a(T) VA(T) 



^lo- 



rn 



Q^{T) ^uj^ + iuj-fo{T) - ujI (29) 

The phason damping 70 from qp backscattering is defined by equation (60) in 
[15]. We estimate its value by introducing a cut-off Uc <t; ksT near the band 
gap, assuming A > ksT and V}, — v for the qp scattering rates (cf. [15]). This 
gives 



lo{T) 



4A keT 



In 



knT 



exp( 



knT' 



(30) 



Thus r = -fo/ujLo = ciu\l-a{iyJkBT)\exp{-A/kBT)/{vFkBTVG) < 1 holds and 
direct phason damping is not important in the frequency and temperature range 
under study. The phason-polariton damping stems from qp dissipation expressed 
by uji. 

The temperature dependence of eA(T) cannot be neglected even when A = Aq 
is considered. It is contained (cf. [15]) in 



eA(T) 



Ao 



/ Ao \' 



The quantity z is defined by 



9i{T). (31) 
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For a qualitative study of finite temperature effects we adopt a simple model in 
which A(T) is fixed at Aq. This is not unreahstic for CDW up to 20% qp fraction 
{N = 0.8). As in [15] we then find 



m 



COLO 



3/2 



(Tat 



rp . 3/2 



N)/gi. 
(33) 

In (33), (Tat denotes the normal state conductivity. For a reference temperature 
of To = Ao/ks = lOOOi^, the parameter Tq is estimated to be of order 10, the 
value we use in the figures. The expansion of the qp fraction 1 — iV for A 3> ksT 



is 



1 - = ^47rV^exp(-2V^), 
and for u — i/b, the backflow parameter bo is approximated by 

1 



bo 



2^/i 



(l-N). 



(34) 



(35) 



5.1 Guided Waves at Finite Temperatures 

Using the same scaling as for T — 0, the matching condition (14) can be written 
for region I (p = 0, 1, 2,3,...) as 



DiT) 



[Z-W) 



a{T) -W _ il{T) 
W -Wo 



71 



+ arctan 



W 



(36) 



^ g{T){hZ - W) [{a{T) - W)/{W - Wo) - il{T)/VW 

Here, the principal values of the square roots and the arctan must be taken. 
The continuation of the p — solution into region 1 gives the surface phason- 
polariton from the former tanh-equation. The following parameters which extend 
zero temperature quantities are used in (36): 



D{T) = (?i(T)Do, (37) 
9{T) = 9i{T)go, 
a{T) = {l-bo{T))yg,{T). 
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The quantity a reduces to unity for T — > 0. It is a measure of the change of 
^Lo{T) with temperature. For small qp concentrations, the decrease of gi{T) 
with increasing temperature makes ujlo{T) larger than the zero temperature 
value lulo- At even higher temperatures, the backfiow parameter 60 but also 
the increased damping 1{T) and the decreased order parameter A(T) become 
relevant for the effective value of u;lo{T): The LO-phason softens and becomes 
overdamped. We do not display this region because it is not properly covered by 
our temperature model. 
The non-retarded version of this equation is 



Zgi 



W 



W-Wo 



TT 



p— + arctan 



hg 



a-W I 
— I 



W -Wq 



(38) 



Corresponding results are shown in Fig. 7. It is seen that the fundamental 
mode {p = 0) is shifted to higher frequencies when the qp fraction increases. 
This reflects the increase of the LO-phason frequency which has been described 
above. The guided waves remain underdamped but the imaginary parts of its 
frequencies do become significant for N > 0.9. The higher order modes (p > 0) 
are much less affected by qp. 

We like to point out that the surface and bulk phasons for c — 00 studied in Sec- 
tion 7.2.3 of [15] appear in a different range than the present phason-polaritons, 
namely at much smaller frequencies, larger wave numbers, and higher tempera- 
tures {ui > ojlo)- Neglect of internal CDW-dispersion is then not permitted. 



5.2 Surface Phason— Polaritons at Finite Temperatures 

Finally, we study the surface phason-polaritons in region 1 starting from the 
appropriate formulae for region 1: 



ir - z r cothX(r) 

^ g{T){hZ - W) [{a{T) - W)/{W - Wq) - il{T)/y/w] ~ \ tanhX(r), 



(39) 



X{T) = D{T) 



\ 



{W -Z) 



a{T)-W _ il{T) 
W -Wq 



Unfortunately, these transcendental equations are difficult to solve, even on a 
computer. Therefore, we limit ourselves to the half space case D{T) — > 00 when 
a simple fifth order polynomial in the frequency u oc \/W remains: 



g(T){hZ - W) [VW(a(T) -W)- il(T){W - Wo)\ = VW(W - Wo)(W - Z). 

(40) 
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Figure 7: Non-retarded limit of polaritonic response of CDW films at finite tem- 
peratures. The wave number scale is Z = e/\q^d'^ /et (2d: film thickness) and W 
means {Reoj /ulo)"^ ■ ^ is the condensate density and guided wave damping is 
displayed as Imuj/ujLo < 0. The quantities go and h are the high frequency CDW 
dielectric constants along (e/\) and transverse to the chains (et) in relation to the 
embedding medium dielectric constant eoo, respectively. 



Its surface phason-polariton solutions are shown in Fig. 8 for several qp densi- 
ties. In contrast to usual surface phonon-polaritons, the qp cause a significant 
temperature variation of the surface polariton frequencies. The surface phason- 
polaritons are underdamped for the qp densities considered. 

6 Summary 

The strong anisotropy of the polar crystal like optical response of quasi one- 
dimensional CDW has a distinguished influence on the coupled excitations of the 
electromagnetic field and CDW phasons. These excitations are grouped into sur- 
face phason-polaritons and guided wave modes and their frequencies are consider- 
ably smaller than the corresponding excitations involving optical phonons. This is 
due to the large effective mass of CDW which leads to a low frequency of the lon- 
gitudinal optical phason. At finite temperatures, quasi-particles cause damping 
and considerable level shifts of both surface phason polaritons and guided wave 
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Figure 8: Surface phason-polaritons in a thick CD W film at finite temperatures. 

W means {Reu /ojloY ■ ^ ^-^ condensate density and polariton damping is 
displayed as bimuo / uolq < 0. The quantities Qq and h are the high frequency 
CDW dielectric constants along and transverse to the chains in relation to the 
embedding medium dielectric constant, respectively. 



modes. We have given a complete discussion of these surface phason-polaritons 
and guided wave modes for the case of P-polarization and for the conducting 
CDW axis lying inside the film. 
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